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ABSTRACT 



A new class of disk MHD equilibrium solutions is described, which is valid within 
the standard local ("shearing sheet") approximation scheme. These solutions have the 
following remarkable property: velocity streamlines and magnetic lines of force rotate 
rigidly, even in the presence of differential rotation. This situation comes about because 
the Lorentz forces acting upon modified epicycles compel fluid elements to follow magnetic 
lines of force. Field line (and streamline) configurations may be elliptical or hyperbolic, 
prograde or retrograde. These structures have previously known hydro dynamical analogs: 
the "planet" solutions described by Goodman, Narayan, & Goldreich. 

The primary focus of this investigation is configurations in the disk plane. A related 
family of solutions lying in a vertical plane is briefly discussed; other families of solutions 
may exist. Whether these MHD structures are stable is not yet known, but could readily 
be determined by three-dimensional simulations. If stable or quasi-stable, these simple 
structures may find important applications in both accretion and galactic disks. 
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1. INTRODUCTION 



The essential role of magnetohydro dynamics (MHD) in understanding the behavior 
of accretion disks is now widely recognized (e.g. Papaloizou & Lin 1995, Balbus & 
Hawley 1998). Usually, however, detailed solutions to the dynamical equations are quite 
complicated: one must contend with turbulence in weak field systems, and exact equi- 
librium solutions, when they can be found at all, are highly technical and rather special. 
Generally, differential rotation will complicate matters by producing a time-dependent 
toroidal magnetic field when a radial field is present. 

Given the above, it might appear that any true equilibrium MHD configuration would 
require at the very least a vanishing radial field component. This is not so. In this 
Letter, we present exact, two-dimensional solutions to the MHD equations, using only 
the approximation of the standard shearing sheet model. (In essence, curvature effects 
are neglected). The solutions are noteworthy for the fact that the field lines, which do 
indeed have radial components, rotate rigidly, even if the disk is not rotating uniformly. 
Furthermore, the field configurations have very simple geometries: hyperbolae and ellipses. 
The dynamical and induction equations are mutually consistent because magnetic forces 
cause departures from the standard epicyclic paths, departures which in these solutions 
lead fluid elements directly along the field lines. Gas pressure gradients are unimportant. 

Of dynamical interest on their own, these solutions may also be astrophysically rel- 
evant. Since only the shearing sheet approximation is used, both galactic as well as 
Keplerian disks are possible venues. The crucial question is whether the solutions are 
stable-or more generously, whether they are stable enough. We do not attempt to answer 
the stability question in this paper, which can be best addressed by three-dimensional 
numerical simulations. It should be noted, however, that the configurations are not 
obviously unstable to the weak field instability discussed by Balbus & Hawley (1991), 
since the field strengths involved fall at the edge of the stability domain. If they are 
long lived, the elliptical solutions (which can be both prograde and retrograde), would 
represent coherent disk structures. One might then speculate that in galactic disks these 
magnetically penned regions in galactic disks become natural sites for molecular cloud 
complex formation. But whether nature makes use of these solutions or not is at present 
unknown. 

Although very different in their detailed physics, these MHD solutions are kin to 
the shearing sheet planet solutions discovered by numerically by Hawley (1987), and 
elucidated analytically by Goodman, Narayan, & Goldreich (1987, hereafter GNG). (The 
GNG solutions allowed for the presence of significant pressure gradients in the underlying 
equilibrium disk; the appropriate limit for comparison with this work is restricted to 
Keplerian rotation profiles.) To understand the similarity, note that the gravitational 
tidal force varies linearly with radius (to lowest retained order). This drives a velocity 
response also linear in spatial coordinates. The associated streamlines are in both cases 
the familiar conic sections. What is surprising, however, is that the presence of Lorentz 
forces can be so easily accommodated by this remarkably simple scaling. 

In § 2, we present a detailed derivation of the MHD structures, and in § 3 a brief 
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discussion of the astrophysical implications of these solutions is presented. 



2. COHERENT MHD STRUCTURES 

2.1 Basic Equations 

Consider a differentially rotating disk, with a cylindrical coordinate system (R, <fi, z) 
centered on the origin. The disk's angular velocity is given by Q(R). The fundamental 
MHD equations are mass conservation 

d\np „ . 

+ V-w = 0, (2.1) 



dt 



the dynamical equation, 



^-K p+ S-^ +( S- v)B ' (2 - 2) 



and the induction equation, 



^=Vx(wxB) (2.3) 



Here, d/dt is a Lagrangian time derivative: 

d d _ 
Jt = d-t +W ' V > 

w is the velocity in the inertial frame, P represents the gas pressure, $ is the central 
gravitational potential, and the other symbols have their standard meanings. We shall not 
require an internal energy equation in our analysis. 

We work in the local (or 'shearing sheet') approximation. This consists of choosing a 
fiducial radius Rq rotating at angular velocity fio? an d erecting a local Cartesian system: 

£ = #-#0, y = Ro(<f>-Clot), z = z. (2.4) 

In the local approximation, x <C Rq, and v x ,v y <C Ro&o, where the V{ are velocities relative 
to a corotating origin. The Alfven velocity is defined by 



v A 



We shall assume initially assume that the Alfven velocity components va x and va v are 
-Ro^o? a result which will prove to be self-consistent. When present, v z and va z will 
also be assumed to satisfy this inequality. 

Let the global scale angular velocity given by a power law of the form Q(R) = 
Qo(Ro/R) q . Keplerian flow corresponds to q = 3/2, a flat galactic rotation curve to 
q = 1. Effecting our coordinate transformation (2.4) and letting v now represent the 
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velocity relative to our new axes leads to the Hill equations (e.g. GNG; Hawley, Gammie, 
& Balbus 1995): 

^+V-v = 0, (2.5a) 
dV - -Iv (P + — ) + _ 20xv + 2qtfxe x - tfze z , 



(it /9 \ 8-7T / 47rp 

f =Vx(vxB), (2.5c) 

where is now measured relative to RQq, and all differential operators may be taken to 
be Cartesian. We have dropped the subscript from O m equation (2.5b). 

2.2 Solutions 

As our first case, let p and P be functions of z only, which is certainly consistent with 
the local approximation. We seek solutions to equation (2.5) in which divergence-free fluid 
velocities and the magnetic field are confined to the xy plane, and va depends upon x and 
y only. In component form, the steady-state mass and induction equations may be written 

diVi = 0, (2.6a) 

VidiVAj = VAidiVj, (2.66) 

with i and j taking on the values x and y. The summation convention on repeated indices 
is used unless otherwise stated. The dynamical equations become 



VidiV x - 2Qv y = -d x j + VAidiVAx + 2qO x, (2.6c) 

VidiV y + 2Vtv x = -d y (^Y^j + v AidiVAy, (2.6d) 
dz (H+^-)=0. (2.6e) 



where the enthalpy H is defined by dH = dP/p. The notation <9j denotes the partial 
derivative with respect to the ith Cartesian coordinate. 

The z equation may be immediately integrated, 

H-H c =^ (2.7) 

where H c is the central midplane value of the enthalpy. The vertical structure decouples 
completely from the planar dynamics of this problem, and we need not pursue it further. 

2.2.1 Hydrodynamic Limit. 
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To orient ourselves, let us recover some familiar solutions for nonmagnetized disks from 
the system (2.6). Setting v x = Ay, v y = Bx (A and B are constants to be determined), 
and va = 0, equations (2.6c,d) lead immediately to 



B(A - 20) = 2g0 2 , A(B + 20) = (2.8) 

There are two distinct branches of solutions. One possibility is A = 0, B = — qQ, which 
corresponds to simple differential rotation of the background flow. The other possibility is 
B = —20, A = 0(2 — q). This is epicyclic motion about the guiding center R = R . The 
streamlines are retrograde ellipses with minor to major axis ratio yjl — q/2. We emphasize 
the simple but important point that fluid elements in the neighborhood of the corotation 
point do not separate from one another as time goes on, but orbit instead about the local 
origin. Therefore, in the absence of Lorentz forces, magnetic field lines placed along the 
epicyclic ellipses would remain undisturbed. 

When q = 2, the ellipses become degenerate, and shrink to straight lines. If q > 2, 
this streamlines become hyperbolic, and correspond to Rayleigh unstable flow (angular 
momentum decreases outward). 



2.2.2 Elliptical Configurations. 

When magnetic fields and Lorentz forces are self-consistently included, it is a remark- 
able fact that equilibrium solutions are still possible. Equations (2.6a-e) have the following 
exact solution: 

v x = Ay, v y = Bx, v Ax = ay, v Ay = (3x, (2.9) 
where the constants A, B, a, (3 are given by 



B = ±Q\I 7-^, A/B = a/ (3 = -e 2 (2.10a) 



and 



Its physical significance of e < 1 becomes evident from consideration of the equations for 
the flow streamlines which coincide with the magnetic field lines: 



x 2 

y 2 + — = constant. (2.11) 

Thus, e is the minor to major axis ratio of the elliptical streamlines. Unlike the hydrody- 
namical case, its value is not fixed by the background Q(R). Furthermore, both retrograde 
and prograde motion is possible; the former (latter) corresponds to taking the — (+) sign 
in equations (2.10a,b). The sign of a has no dynamical consequence (it determines the 
sense of the current flow), but once chosen, (3 must be of the opposite sign. Note that 
that the magnetic field lines do not shear, even though a radial component is present and 



6 



there is background (global) differential rotation. Rather, the field lines precisely track the 
velocity streamlines, and rigidly rotate with the modified epicycle. The "miracle" is that 
this behavior is dynamically fully self-consistent. 

Let us next consider the conditions under which retrograde or prograde ellipses are 
obtained. Retrograde motion corresponds to A > 0, B < 0, i.e., to taking the minus sign 
in equations (2.10a,b). The requirement that a 2 > means that 

1 - q/2 < e 2 < 1 (retrograde). (2.12) 

For a Keplerian disk, e must therefore exceed 1/2, for a galactic q = 1 law, e must 
exceed \[2/2. The Keplerian planet solutions of GNG, which were also retrograde ellipses, 
had a complementary restriction on the domain e: < e < 1/2. What is inaccessible 
hydro dynamically becomes available magneto-hydrodynamically, and vice-versa. Note 
that when 1 — e 2 approaches q/2, we recover the the hydrodynamical solution of the 
previous subsection. This is what ultimately limits the domain of e for retrograde ellipses. 
The inclusion of magnetic hoop stresses "plumps" retrograde ellipses into more circular 
structures. 

Prograde motion corresponds to A < 0, B > in equation (2.7). The prograde branch 
is not restricted by the condition a 2 > 0, since this is obviously guaranteed by choosing 
the plus sign in equation (2.10b). Hence, the domain of e for prograde ellipses, 

< e 2 < 1 (prograde). (2.13) 

extends beyond the retrograde domain. 

The vorticity of the ellipses is given by 



B - A + 20 = O 



2±(l + e 2 ) 



2q 



1 -e 



Equation (2.12) implies that retrograde ellipses always have negative vorticity; prograde 
ellipses have positive vorticity. 

The Lorentz forces are — (3 2 (1 + e 2 )x in the x direction, and — e 2 /3 2 (l + e 2 )y in the y 
direction. This is an inwardly directed force (confining) for both prograde and retrograde 
ellipses, for a given e larger in magnitude for the prograde solutions. Prograde rotation 
has greater associated vorticity than retrograde rotation, and larger confining forces are 
required. 

2.2.3 Hyperbolic Configurations. 

Hyperbolic configurations are another possible local solution in a disk. Following the 
ABa(3 parameterization of equation (2.9), we obtain 



S = ± ^TT1^' A I B = *I$ = M ^ ( 2 - 14 ) 
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and 

1-M 4 \ \j q/2 J 

As before, the sign of a is not dynamically significant, but once chosen, a and (3 must have 
the same sign. These solutions have flow streamlines and magnetic field configurations 
given by the hyperbolae 

x 2 

y 2 — = constant. (2-15) 

Prograde hyperbolic flow corresponds to B > 0, M 2 < 1; retrograde flow to B < 0, 
M 2 > 1. (Note that with a 2 is always positive under these restrictions.) Lorentz forces 
for hyperbolic flow are directed neither radially inward nor outward (relative to our local 
origin), but vary with quadrant. 

2.2.4 Constant Density Solutions. 

The planar field configurations of the previous section are fully compatible with 
vertical stratification in the disk. If one further restricts the disk structure to the case of p 
and P being independent of z, a greater variety of possible field configurations exists. This 
is not an uninteresting limit, because numerical simulations are often performed ignoring 
the vertical disk structure, and because the local midplane structure of astrophysical disks 
is approximately one of constant density. (More precisely, p = po — 0(z 2 ), and we shall 
work to 0{z).) 

We may first note that if p and P assumed to be constant, nothing in our original 
solution would change if a constant vertical magnetic field were present. This is one possible 
generalization. But it is also possible to find qualitatively new solutions. 

If z components of the velocity and magnetic field are present, but vertical structure 
is absent, instead of equation (2.6e), our new z equation is 

VidiV z = -d z {v\/2) + VAidiVAz (2.16) 

Consider now a velocity of the form 

v x = Az + D, v y = Bx, v z = Cx, (2-17) 

where A, B, C, and D are constants, and an Alfven velocity of the form 

v A = 7v (2.18) 

with 7 constant. The equations (2.6a,b) are automatically satisfied, while equations 
(2.6c,d) and (2.16) lead to elliptical and hyberbolic streamline solutions, as before. El- 
liptical solutions take the form 

C 2 = ^ 2 [2(1 + e 2 ) 2 - ? ], C/A = -e 2 , 7 2 = ^ , (2.19a) 
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and 



B = -20(1 + e 2 ). 



(2.196) 



These solutions have the interesting property that D is completely unconstrained, and the 
entire ellipse can stream, in bulk, in the radial direction! Unusual streaming behavior 
is also characteristic of the magnetorotational streaming solutions in two-dimensional 
axisymmetric simulations (Hawley & Balbus 1992). However, these streaming solutions are 
known to be unstable to Kelvin-Helmholtz types of instabilities when three-dimensional 
structure is permitted (Goodman & Xu 1994, Hawley, Gammie, & Balbus 1995). The 
ordered flow breaks down into MHD turbulence. 

3. DISCUSSION 

The solutions presented in this paper represent both O-type (elliptical) and X-type 
(hyperbolic) neutral points in the local magnetic field topology. In essence, we have shown 
that each of these topologies has a range of well-defined field strengths that leads to exact 
static equilibrium solutions to the local MHD equations in a shearing disk. It is remarkable 
that such solutions exist in the presence of differential rotation. 

Unless e 2 — 1 is very small, the local Alfven velocities of the solutions are of order 
rO, where r is the radial distance from the local origin. If r < the disk scale height, the 
magnetic field strengths are comparable to or less than thermal values. As e — > 1, field 
strengths becomes unbounded, and the configurations are probably unstable. At more 
nominal field strengths, the question of stability is unclear. Both the Parker (1966) and 
magnetorotational (Balbus & Hawley 1991) instabilities are potentially disruptive, but in 
either case the lengthscale and magnitude of the magnetic field make direct application of 
the relevant criterion marginal. The best option at this point is direct numerical simulation. 

But these concerns should not stifle speculation. Even if these structures prove 
unstable on rotational time scales, they may still be astrophysically interesting. X-type 
neutral points facilitate reconnection, and their presence as an equilibrium configuration 
at thermal field strengths may be an important saturation mechanism for accretion disk 
dynamos. Another application of interest is the interface between an accretion disk and 
the magnetosphere of the central source. The magnetic field enforces rigid rotation with 
a fiducial radius, the corotation radius, while differential rotation is maintained outside 
corotation. This work opens up the possibility that magnetic fields at the interface need 
not be sheared! The open field lines of our hyperbolic solutions may globally join the 
magnetospheric fields to the disk interior without producing large off-diagonal Maxwell 
stresses. 

O-type elliptical structures offer interesting possibilities for galactic disks. If we 
identify these large scale vortices with sites of molecular cloud complexes, there are a 
number of straightforward predictions. We should see, of course, at least a vaguely elliptical 
morphology, with the major axis oriented azimuthally, and magnetic fields following this 
morphology. Both retrograde and prograde rotation of the complexes is possible. The 
minor-to-major axis ratio of retrograde systems must exceed 0.7; however prograde systems 
with considerably smaller values should exist. 



9 



Why have numerical simulations done to date not shown these structures? Three- 
dimensional MHD disk simulations have yet to reveal the midplane elliptical and hyperbolic 
configurations. There are several possible reasons for this. One, unlike the GNG solutions, 
the density does not fall to zero at the edge of an ellipse (or at some convenient point 
in an hyperbola). This means these solutions must be actively confined by the ambient 
disk, and conditions for this need not arise spontaneously. Two, numerical simulations 
have concentrated on subthermal magnetic fields, because their presence leads to disk 
turbulence. If the field saturates at subthermal values, and it generally does, the O and 
X neutral points may not be able to form. (In the simulation where the traveling ellipses 
were seen, the field strength grew to suprathermal values.) The best strategy might be 
to start with thermal fields and allow some combination of differential rotation, internal 
dynamical instabilities, and external driving to a new field configuration. Finally, there is 
always the possibility that they structures never form spontaneously because they are too 
unstable. 

Undue pessimism, however, is not yet warranted. Properly crafted three-dimensional 
simulations should be able to clarify the most important uncertainties. 

Part of this work was completed when one of us (SAB) was a visitor with the radio 
astronomy group at the Ecole Normale Superieure, and he would like to thank E. Falgarone 
and M. Perault for their generous hospitality and advice. This work has been supported 
by NASA grants NAG5-4600, NAG5-3058, and by NSF grant AST-9423187. 
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